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Abstract. The gravitational recoil or "kick" of a black hole formed from the merger 
of two orbiting black holes, and caused by the anisotropic emission of gravitational 
radiation, is an astrophy sic ally important phenomenon. We combine (i) an earlier 
calculation, using post-Newtonian theory, of the kick velocity accumulated up to the 
merger of two non-spinning black holes, (ii) a "close-limit approximation" calculation of 
the radiation emitted during the ringdown phase, and based on a solution of the Regge- 
Wheeler and Zcrilli equations using initial data accurate to second post-Newtonian 
order. We prove that ringdown radiation produces a significant "anti-kick" . Adding 
the contributions due to inspiral, merger and ringdown phases, our results for the net 
kick velocity agree with those from numerical relativity to 10-15 percent over a wide 
range of mass ratios, with a maximum velocity of 180 km/s at a mass ratio of 0.38. 
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1. Introduction and summary 

The gravitational recoil of an isolated system in response to the anisotropic emission 
of gravitational radiation (sometimes also called the "kick") is a phenomenon with 
potentially important astrophysical consequences [T]. One of the most intriguing is 
the possibility that a massive black hole formed from the inspiral and merger of two 
progenitor black holes could receive enough of a kick to displace it from the center of 
the galaxy where the merger occurred, or to eject it entirely from the galaxy. This could 
affect the growth history of massive black holes [2j. Observational evidence for such a 
kicked black hole has even been reported [3]. 

The calculation of such kicks within general relativity has been carried out in a 
variety of ways. Earlier analytic or semi-analytic estimates of the gravitational recoil 
include a perturbation calculation (valid for small mass ratios) during the final plunge 
[I], a post-Newtonian calculation valid during the inspiraling phase together with a 
treatment of the plunge phase [5], an application of the effective-one-body formalism 
[6], and a close-limit calculation with Bowen-York type initial conditions [7j. 

Following recent advances in numerical calculations of binary black holes [HI El EI] , 
the problem of gravitational recoil received considerable attention from the numerical 
relativity community. These computations led to increasingly accurate estimates of the 
kick velocity from the merger along quasicircular orbits of black holes without spin 
PH CGI H3 M, ES], and with spin HQ HH HBJ; from head-on collisons [Tg]; and from 
hyperbolic orbits [20]. In particular these numerical simulations showed that very large 
kick velocities can be obtained in the case of spinning black holes for particular spin 
configurations. Nevertheless, as the very detailed multipolar analysis of the binary black 
hole recoil by Schnittman et al. [21] illustrates, analytic and/or semi-analytic methods 
are still very useful for gaining more physical understanding of the relaxation of binary 
black holes to their final equilibrium state. 

In the simplest case of unequal mass, non-spinning black hole binaries on 
quasicircular orbits, the kick velocity as a function of time shows a very distinctive 
pattern [121 E3! : the recoil increases monotonically during the inspiral and plunge 
phases up to a maximum around the onset of merger, and then decreases quickly to 
a final asymptotic value, as much as 30 percent smaller than the maximum. This 
braking occurs during the phase where the newly formed black hole emits gravitational 
radiation in a superposition of quasinormal "ringdown" modes, and is known as the 
"anti-kick". For a reduced mass parameter value r\ = m\m<il{vii\ + m2) 2 ~ 0.19, the 
value for which the kick is a maximum, the peak value is around 250 km/s while the 
final kick is around 175 km/s. 

Building on previous work based on a multipolar post-Minkowskian formalism 
[221 1231 121!, Blanchet, Qusailah and Will [5] (hereafter BQW) derived the linear 
momentum flux from compact binaries at second post-Newtonian (2PN) order beyond 
the leading effect. BQW augmented their 2PN estimate of the recoil up to the innermost 
circular orbit (ICO) by integrating the resulting 2PN-accurate flux from the ICO down 
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Figure 1. Left panel: Comparison with numerical calculations jTTJ [T^l EH HH ITS] 
and other analytic or semi-analytic methods [5J El E] • Right panel: Comparison with 
fitting formulas derived from this paper (black) and from numerical relativity results 

to the horizon on a plunge geodesic of the Schwarzschild geometry. They found that 
the recoil monotonically increases as the plunge progresses. Within their error bars, the 
resulting recoil was found to agree well with the maximum value of the kick velocity 
as calculated by numerical relativity up to the onset of the anti-kick. And indeed the 
maximal kick velocity in numerical computations occurs more or less at a separation of 
roughly 2M, where M = m,\ + m 2 is the total mass, as inferred from the times at which 
the maximum kicks were found to occur in [T2l WQ (such an inference cannot be made 
rigorously, of course, but does coincide roughly with where other numerical diagnostics 
indicated the onset of the merger). At this point, the BQW computation ended for lack 
of a method to evaluate the contribution from the subsequent ringdown phase. This 
paper reports the results of incorporating such a method. 

Le Tiec and Blanchet [27] have developed a "close-limit approximation" (CLA) 
for black hole binaries that uses 2PN-accurate initial conditions. In this 2PN-CLA 
framework, the 2PN metric for two bodies in close proximity is recast as a perturbation 
of a Schwarzschild black hole. The resulting perturbation is then used as initial data 
to evolve numerically the Regge- Wheeler and Zerilli equations in order to calculate the 
gravitational radiation emitted subsequently. The purpose of the present paper is to 
use the resulting waveforms to compute the recoil generated during the ringdown phase. 
By adding vectorially the 2PN-CLA results to the 2PN results of BQW for the inspiral 
plus merger, we prove that the effect of the ringdown on the recoil is indeed to produce 
an anti-kick, and we find that the total kick generated by the inspiral plus merger plus 
ringdown phases is in good agreement with numerical computations for non-spinning 
black hole binaries. 

Our central results are shown in Fig. [T] In the left panel, the curve with error 
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bars (green in the color version) shows our combined kick velocities, with the error bars 
estimated by varying the radius at which the 2PN and CLA methods are matched. The 
top curve (red) is the kick from the BQW pure 2PN calculation up to the merger. The 
sequence of dots and accompanying dashed lines (blue) are from an exhaustive series of 
numerical simulations by Gonzalez et al. [Ml [T5] . In this figure are also shown individual 
points and error estimates from some earlier analytic or semi-analytic estimates such as 
[HI E], as well as other numerical computations. 

A number of authors have fit the kick velocity to the empirical formula p2, E] 

V &t = Ar] 2 (l - Ar]) 1/2 (1 + B rj) x 10 3 km/s. (1) 

The leading rf{\ — Arj) 1 / 2 dependence derives from the lowest-order, or "Newtonian" 
calculation (281 EE]. A fit to our 2PN-CLA results yields A = 9.5 and B = 0.3. 
The right panel of Fig. [T] shows the ratio of our kick velocities to this fitting formula 
(black), together with similar ratios to fitting formulas derived from numerical relativity 
[251 126] . It can be seen that our kick velocities are systematically higher than those from 
numerical relativity in the equal-mass limit, and systematically lower in the small mass- 
ratio limit. However, we find it striking that in the regime where the kick velocity 
is substantial, say between r\ = 0.08 and r\ = 0.24, our 2PN-CLA calculation agrees 
with numerical relativity to 10-15 percent, and confirms the intuition that ringdown 
radiation generates an anti-kick that partially offsets the kick accumulated during the 
plunge. 

The remainder of this paper is organized as follows: In Sec. [2] we describe briefly the 
2PN and CLA methods for calculating the gravitational recoil from the inspiral, plunge 
and ringdown phases (full details will be presented elsewhere |30j), and we describe 
the numerical implementation of the method and some of the checks and diagnostics 
performed. Concluding remarks are made in Sec.[3j We use geometrical units G = c = 1. 

2. Gravitational recoil from the inspiral, plunge and ringdown phases 

The flux of linear momentum carried away by gravitational waves from a general isolated 
system can be written in terms of the gravitational- wave polarization states h + and h x 
as [31] 

^= lim (fn i \h+-ih x \ 2 dn\, (2) 

dt r-»+oo [ lb7T J J 

where the dot stands for a partial time derivative, and dfl is the solid angle associated 
with the direction of propagation n l . 

BQW [5j expanded the waveforms in terms of radiative multipole moments using the 
post-Minkowski framework, and expressed them in terms of source multipole moments 
of mass-type and current-type, valid to 2PN order, including contributions of "tails". 
Restricting to binary systems on quasicircular inspiral orbits, they obtained dP l /dt, and 
after integrating with respect to time, dividing by the total mass M and changing sign, 
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obtained the recoil velocity 
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where x = (Ma;) 2 / 3 , uj being the orbital angular frequency, and u % is the unit 
vector directed from the less massive toward the more massive body (See [32J for 
a generalization including spin effects.) This formula gives the kick velocity V™ spiral 
accumulated during inspiral up to the ICO (defined by xico — §)■ 

Starting with V^ spiral , which is always a small contribution to the total kick, BQW 
then integrated the 2PN expression for dP*/dt along a plunge orbit of a Schwarzschild 
black hole of mass M from the ICO down to a radius of order 2M. A key to that step 
was to change integration variable from coordinate time t, which is singular on the event 
horizon, to a "proper angular frequency" variable to = dtp/dr of the plunge orbit, which 
is regular on the horizon. The resulting net kick was the data for the curve labeled 
BQW (red) plotted in Fig. [I] Here we repeat this calculation, except that we terminate 
the plunge integration at a Schwarzschild coordinate radius r matc h, whose value is chosen 
to lie between 2M (the minimum allowed by the method) and 2.5M. At this radius we 
match the 2PN kick (namely V^ spiral + Vp\ unge ), to the result of our 2PN-CLA calculation, 
to which we now turn. Later we will test the sensitivity of the final result to the value 
of the matching radius r match . 

The idea of the 2PN-CLA method [27] is to take the spacetime metric for a binary 
system accurate to 2PN order, where the two bodies are on a quasicircular orbit of initial 
separation r\2 (in Schwarzschild-like coordinates), which is assumed to be of order 2M. 
The metric is then re-expanded in powers of r 12 , resulting in a Schwarzschild metric of a 
black hole of mass M plus correction terms that vanish in the limit 77 — > 0. Carrying out 
a multipolar expansion, one can identify the (£, m) components of the Regge- Wheeler 
and Zerilli functions used in black-hole perturbation theory. 

It is possible to express the linear momentum flux (J2J) in terms of the Regge- Wheeler 
and Zerilli functions as defined in Eqs. (5.1) of [27J. The waveform takes the form 

K-iK = l -Y. (*& + i *S) -*nm + 0{r->) , (4) 

where the superscripts (e) and (o) denote even and odd-parity respectively, the 
summations on the integers £, m range from 2 to infinity for £, and from — i to i for m, 
and where -2Ye tm (9, ip) are the spin-weighted spherical harmonics of spin —2 [33J, 134] . 
More details on this standard result are given in the Appendix of [27] , including its less 
well-known generalization at any order in r~ l . We insert Eq. Q into Eq. ^ and find 
(see [35] for details) 

dP x .dP, 



dt dt 8?r „ 



. dfy 1 • T fe) t W , ( T (g) T T (O) t (°) 



(5) 
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Figure 2. Magnitude of the recoil velocity Vri ng down generated during the ringdown 
phase as a function of the symmetric mass ratio rj for different initial separations ri2 . 



(£+3)! / (l+m+l)(l+m+2) 



where a^ m = 2(£ - 1) (£ + 2) ^(1 - m)(^ + m + 1) and =■ (m)(<? _ 2)! y (2*+i)(2M-3) 
and the overbar denotes complex conjugation. Because of the symmetry with respect 
to the orbital plane, we naturally find dP 2 /dt = 0. This would no longer remain true if 
we were to include spin-orbit coupling terms for spinning black holes in the initial PN 
metric. The master functions obey the wave equations 

(#-£ + V<^>) *<*? = (), (6) 

where the tortoise coordinate r* is related to the Schwarzschild radial coordinate r by 
r* = r + 2Mln(r/2M — 1), and where the potentials V^ 6 ' '* are given by 

^=(l-^)(^-^>), (7) 

with Uf = MA ';yS;r M) and W<°> = l, where A, = \(t - + 2). 

We evolve the equations ^ with 2PN-accurate initial conditions computed with 
an initial separation ri 2 , and developed in the CLA as described above; detailed inputs 
are Eqs. (3.6)-(3.7) and (4.4)-(4.5) of Ref. [27] (see also [7J for an alternative CLA 
calculation using different initial conditions). Then, inserting the numerically generated 
master functions ^q- (5)> we calculate the flux of linear momentum up to 

octupolar order during the ringdown phase, and, integrating that with respect to time, 
dividing by M and changing the sign, we obtain the ringdown contribution V r \ ngdown to 
the total kick. Fig. [2] shows the magnitude K. ing down = I Kingdown I as a function of 77 for 
various values of 7*12. 

Choosing r 12 = r match , we then add up vectorially the results for the inspiral, plunge 
and ringdown phases to obtain V % = V^ spiral + V^\ unge + V7 ingdown . In all cases we find that 
the direction of the ringdown kick is approximately opposite to that of the accumulated 
inspiral plus plunge kick. Not surprisingly, the direction or phase of the inspiral plus 
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Figure 3. Dependence of the result for the direction (left panel) and magnitude (right 
panel) of the recoil velocity upon the matching radius r matc h. 



plunge kick is sensitive to the radius r matc h at which the 2PN calculation terminates. 
But, most satisfactorily, when we add the ringdown kick, the final direction is relatively 
insensitive to the value of r matc h, as shown in the left panel of Fig. [3] (see especially 
the inset panel). Similarly, the right panel shows that the magnitude of the total kick 
velocity is also only weakly sensitive to r matc h. In Fig. [TJ we have chosen r matc h = 2.2M 
as being a value where the phase and magnitude of the kick do not vary too much, and 
estimated error bars by varying r matc h between 2M and 2.5M. Unlike BQW, we have 
not attempted to estimate errors caused by the neglect of higher PN corrections in the 
CLA method. 



3. Conclusions 



We have found that the recoil velocity of coalescing, non-spinning black holes can be 
calculated using a combination of post-Newtonian theory for the inspiral and plunge, 
and a close-limit approximation for the ringdown, with results that agree closely with 
those from full-scale numerical relativity. We have also used this method to determine 
the total energy and angular momentum radiated during inspiral, plunge and ringdown; 
details will be published elsewhere [30] . An obvious, though non-trivial next step would 
be to incorporate the effects of spins in this approach. 
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